Phase rigidity breaking in open Aharonov-Bohm ring coupled to a cantilever 
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The conductance and the transmittance phase shifts of a two-terminal Aharonov-Bohm (AB) 
ring are analyzed in the presence of mechanical displacements due to coupling to an external can- 
tilever. We show that phase rigidity is broken, even in the linear response regime, by means of 
\ inelastic scattering due to phonons. Our device provides a way of observing continuous variation of 

the transmission phase through a two-terminal nano-electro-mechanical system (NEMS). We also 

■ propose measurements of phase shifts as a way to determine the strength of the electron-phonon 
coupling in NEMS. 
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Nanoelectromechanical systems (NEMS) have been a subject of extensive research in recent years. The possibility 
I/"") \ of combining electrical and mechanical degrees of freedom on the nanoscale may give rise to technological advantages 
t-H . as well as manifestations of fundamental physical phenomena. From a technological point of view the interest is 
largely due to the many applications that may be realized using NEMS.- It follows ubiquitously from the concrete 
possibility of downsizing devices from micromechanical into nanoelectromechanical systems. Among the many NEMS 
phenomena of considerable physical interest, we focus in this paper on the effect of quantum-coherent displacements in 
the presence of a Aharonov-Bohm effect in a one-dimensional ring symmetrically connected to two external leads. On 
the nanoscale level the mechanical forces controlling the structure of the system are of the same order of magnitude 
as the capacitive electrostatic forces governed by charge distributions. This circumstance is of the utmost importance 
when analyzing electron tunneling. For example, some predictions on the role of electron-phonon interactions on 
the electron tunneling have regarded the modifications of the peak-to- valley current ratio due to inelastic processes 
with important consequences in device applications 2 ^. Another important issue in NEMS physics is related to the 
way mechanical displacements perturb phase coherent charge transport through closed loops as a consequence of the 
shifting of electron trajectories. For example, important information come from the Aharonov-Bohm (AB) oscillations 
in mesoscopic rings^i 6 -! 7 -^. Here, the loss of coherence, or dephasing, due to inelastic scattering caused by phonons 
Ch \ reveals in the suppression oihe/e oscillations. Another important aspect unnoticed in previous studies of NEMS is that 
inelastic scattering breaks phase-rigidity, i.e. the linear conductance is not symmetric in the AB phase. Here we show 
11 ' that the observed asymmetry can be tuned continuously by changing the electron-phonon coupling, demonstrating 
that the phase of the linear conductance in a two-terminal AB interferometer is not rigid when tunneling is assisted 
by phonons. The problem of phase rigidity and its breaking in ring and ring-dot systems has been largely investigated 
both theoreticall y 7 ' 9 ! 10 ' 11 ! 12 ' 13 ! 15 and experimentall y 16 ' 17 ! 18 ' 19 . Here we address for the first time the problem in the 

■ framework of nanoelectromechanical systems and propose a way to detect phase rigidity breaking. 

^v^j ! We consider a one-dimensional ring symmetrically coupled to two leads and to a mechanical cantilever. A coupling 

■ between the electrons travelling in the lower arm of the ring and the cantilever, whose tip is suspended over the arm, 
can be set up by developing a uniform electric field between the tip, the cantilever and the lower arm. Electrons 
couple approximately linearly to the cantilever position, thus leading to a coupling between the flexural phonon 
modes of the cantilever and the local density of the electrons on the ring arm. Furthermore, because the coupling 
strength decays rapidly with increasing frequency, for micron scale cantilevers only the fundamental flexural mode 
is relevant. Therefore, at low enough temperature, the cantilever can be treated as a single quantum mechanical 
oscillator. Of course the dephasing behavior of the electrons due to the cantilever depends on the relative magnitudes 
of the electron dwell time on the lower arm of the ring and the cantilever period. A similar device has been proposed 
by A. MacKinnon and A. D. Armour as a which-path device for electrons [2fj| . 

In the following we generalize the scattering matrix approach in a way suitable for electron dwell times shorter or 
comparable to the cantilever period. In this way both elastic and inelastic contributions to the scattering process will 
appear. The Hamiltonian of the system shown in Fig[T]is written as: 
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H = H r + H C + Hint, (1) 

where Ti. r and Ti. c are the Hamiltonian of the ring (in polar coordinates) and the cantilever, respectively: 



2 



while Hint describes the interaction between the lower arm of the ring (0 < ip < tt) and the cantilever 

H in t(x = R<p,y) = < . 4 

I otherwise, 

where a represents the average of the electron density on the lower arm of the ring times the uniform electric field. 
The interaction between the electrons on the lower arm of the ring and the cantilever is modelled as a linear coupling 
between the displacement of the flexural mode and the average electron density of the lower arm of the ring, since 
when electrons are not in the ring there is no displacement of the cantilever (we assume that all positions are measured 
from the equilibrium height of the cantilever) . The linear coupling is a valid approximation in the limit in which the 
displacement of the cantilever from its equilibrium height is small compared to the characteristic length of the harmonic 
oscillator. The Hamiltonian of left and right lead is given by the free particle Hamiltonian Tli ea d = —h 2 d 2 2 /(2m e ). 
By combining the harmonic potential of the cantilever with the interaction term 7i jri t, the effective potential along 
the lower arm can be rewritten as: 

1 o . ct , o a 2 , . 

V eff = -M c uj 2 (y + - 2 -. 5 

In the following we solve the scattering problem, generalized to include inelastic scattering as in Ref 21. For our 
purposes, we set up a scattering problem by following the method of quantum waveguide transport on network o 22 ' 23 . 
One main problem is the boundary conditions at the intersection with the external leads. In this case the Griffith 
boundary's condition 24 state that (i) the wave function must be continuous and (ii) the current density must be 
conserved. We assume that when an electron moves along the upper arm in the clockwise direction from tp = 0, 
it acquires a phase $ab/2 at the output intersection ip = tt, whereas the electron acquires a phase —$,43/2 in the 
counterclockwise direction along the lower arm when moving from ip = to <p — tt. The wave-function for the right 
(out) /left (in) lead and the upper (up) /lower (low) arm of the ring may be written as: 

00 

ip in (x,y) = (j)i(y)e lkx + r m i<!? m (y)e~ lkmX (6) 

oo 

ipout(x,y) = t »ii ( t ) rn(y)e lkmX (7) 

m=0 

oc 

i>u P (x, y)=J2 4>m(y)(a m ie lk+ ^ x + b ml e~™ f e~ lk ™ x ) (8) 
m=0 
00 

i>lo«,(x, y)=J2 ^(V + yo)(A ml e fk ^ x + B ml e™ f e~ fk ™ x ), (9) 
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where 4> n (y) represents the n— phonon state, and we have defined / = ^ff-- The electron wavevectors, determined by 
the energy conservation, are k rn = yjk 2 + p 2 (l — m), k rn — ^Jk 2 +p 2 (l — m + j 2 ), p 2 = 2m e Lu/h, k^ = k m ± Trf/L, 
&m = i n f I L = being the half circumference) . The phonon wave function in the lower arm is shifted of 



the quantity yo = j^~p and we have defined the dimensionless parameter 7 = yoy The unknown transmission 

and reflection coefficients are determined by solving the following system of equations corresponding to the Griffith's 
boundary conditions: 

ipin (0, y) = ip up (0, y) = ipdow7i(0, y) 

ipaut^R, y) = ipu P (rtR, y) = tpdown^R, y) 

d x tpu P (x = 0,y) + d x ipdown(x = 0, y) = d x ip in (x = 0, y) 

d x tpout(x = ttR, y) - d x tpup(x = ttR, y) = d x ipdown(x = nR, y). (10) 

To eliminate the dependence on y, we project each equations on the phonon state <fis(y)j the only non-trivial projection 
being <p s (y)*(f> m (y + y Q )dy. Since y Q is small quantity we may expand (j> m (y + y ) as <p m (y + y ) rj (f> m (y) + 
"f(y/rn4> m -i(y) — y/m + l(j> m +i(y)), which allows immediately to get the following result: 



4>s(y)*4>m(y + yo)dy w S 8im + j(Vm6 Sim -i - Vm + 15 s ,m+i)- (H) 
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Since 7 is a small quantity, only one-phonon processes are relevant and the scattering problem can be solved with 
arbitrary numerical accuracy. 

In order to solve the linear system (|10| for the unknown transmission and reflection coefficients, a pruning procedure 
has been applied (similar to Refl2lT). The method is based on the simple idea of fixing the number of phonon modes at 
s = S, and than iteratively removing all the coefficients with phonon indices greater than S, adjusting the remaining 
transmission and reflections coefficients in such a way to have a total probability (reflection and transmission) equal 
to one. Once fixed s — S, the problem reduces to solve a linear system in the variables r s i, t s i, a s i, b s i, A s i, B s i, with 
s = 0, . . . , S. In the following we consider S = 4. Once the scattering problem is solved, the transmission probability, 
both for elastic and non-elastic processes, is defined as follows: 

Tmi = \t m i\ 2 -r-, (12) 

where T m i represents the transmission probability from the /-phonon state to the m-phonon state. In what follows 
we assume that the starting phonon state is I = 1, and the initial electron momentum ki is fixed to k. The total 
transmission probability is obtained as ^^i. We solved the scattering problem for s = 0,1,2,3,4 phonons, 
assuming the initial state of the cantilever corresponding to the single phonon state (I = 1). In FigTJl the transmission 
coefficients and the transmittance are shown as a function of the external flux / for kL — 3.5, 7 = 0.2 and pL = 3. 
Even tough the main contribution to the transmission is generally due to the elastic term tn, a strong influence of 
the inelastic term i 2 i is observed close to half-integer values of the flux. The lifting from zero of the conductance close 
to an half-integer flux is the signature of a phonon-assisted tunnelling (PAT). To better discern elastic from inelastic 
contributions to the scattering in FigJ3]the transmission coefficients in, £21,^01 ar e plotted as a function of kL for 
7 = 0.2, pL = 3 and / = 0.49,0.5,0.51 (from top to bottom). As above, a strong competition between elastic and 
inelastic scattering amplitudes is observed close to half-integer values of the flux. A general feature of all the panels in 
Fig|3]is that in the low energy region (kL < 3) the scattering amplitude is dominated by the coefficient toi describing 
the process of emission of a phonon. In the intermediate energy region (3 < kL < 5) the transmission appears 
strongly affected by a resonant peak related to absorption of a phonon (£21)1 while in the high energy region (kL > 5) 
an alternating behavior is observed. Further, in Fig(3Jtop panel), close to kL = 4, only the inelastic coefficient t^i 
contributes to the transmission. This follows from the fact that an electron can be transmitted only by changing the 
cantilever state, i.e. by means of the absorption/emission of a phonon in the final state. Similar behavior can be seen 
for other initial momenta kL. 

Let us apply the general conductance formula in the linear response regime (Landauer-Biittiker formula 25 ), to 
investigate the dependence of the transmission on the flux 4> and through that obtain the behavior of the phase 
shift. The transmission phase through a mesoscopic AB ring has the general property of phase rigidity coming from 
the two-terminal nature of the set-up and which is generally based on time-reversal symmetry-breaking and current 
conservation 2 .-, i.e. the transmission probability amplitude satisfies the property T a p(kL,<f>) — Tp a (kL,—4>) as a 
function of the flux and of the incident electron energy (a, (3 — L, R being the lead index). Combining time-reversal 
symmetry requirement with current conservation \T a p(kL, 4>)\ 2 — \Tp a (kL, </>)| 2 , implies that the linear conductance 
Q is an even function of the flux, whose Fourier transform is 

Q(4>) = Go + ^2G,iCOs(n(f> + 5 n ), (13) 

n 

and 4> = 2tt/. Obviously, the phase shift 6 n can take only two values and 7r— . Thus, the phase of the two 
terminal AB ring has to be rigid, or change abruptly by tt as the accumulated phase in one arm is being varied. This 
peculiar behavior, known as phase rigidity has been studied extensively 7 -! 9 - ' 11 ' 12 ' 13 ' 15 ' 16 ' 19 . Exp er iment all y 1 6 ' 1 7 ' 1 8 and 
theoretically 10 various methods of avoiding phase-rigidity in rings and ring-dots systems were previously discussed. 
Here we are proposing, for the first time, observation of phase rigidity breaking in a nanoelectromechanical system as 
an alternative way of observing the continuous variation of the transmission phase through a two-terminal mesoscopic 
system. Our results in Fig0] show that for different values of the interaction 7 inclusion of inelastic scattering, breaks 
phase rigidit y 28 ' 29 . The first evidence for this comes from the fact that T is non symmetric. This is clearly seen in the 
vicinity of the antiresonances in the conductance, which are approximately located at the energy values determined 
by the electron-phonon coupling 7: 

K l nm L » v /^ 2 (n + /) 2 + (pi) 2 G-m + 7 2 ), (14) 

where n labels the momentum of the eigenstates in the ring, m(l) labels the final(initial) phonon channel. The amount 
of phase rigidity violation depends on strength of 7 and it is related to the momentum of incoming electron. For 
instance, in Figdl the total transmission as a function of the flux / is shown, for pL = 3, kL = 5.5 and varying 7 from 
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zero to 0.2 (from top to bottom with step 0.05). Since the phase shift is a measurable quantity in the experiments, in 
FigsO the behavior of i5i is shown as a function of 7 for different values of the antiresonancc momenta. As shown, in the 
absence of interaction, i.e. under Onsager symmetry, the phase shift is equal to 71", corresponding to a local minimum 
in the conductance versus flux close to / = 0. When the interaction is turned on, the phase shift Si converges quite 
rapidly to the values tt/2 or 37r/2, depending on the value of the momentum kL. The interaction region characterized 
by a rapid variation of the phase shift can be exploited in experiments to obtain information on the electron-phonon 
coupling simply by measuring the phase shift of the linear conductance. Furthermore, since this analysis can be done 
for all the antiresonances, the value of the interaction 7 can be obtained by experimental measurements of the phase 
shift. The phase rigidity breaking we have discussed, has close similarities to that described by Sun et al. in [l5| 
where a two-terminal modified AB ring with a quantum dot inserted in one arm has been considered. There, the 
phase rigidity is broken in the linear regime by applying a time- varying microwave (MW) field on the quantum dot. 
Further, the behavior of the phase shift in FigsOis similar to that shown in Fig. 4 of Ref,19 for non-linear transport 
in an AB ring. Even tough, the experimental analysis as a function of a perpendicular magnetic field is focused on 
nonlinear contribution to the conductance, the analogy in the phase shift behavior can be associated to an effective 
energy shift in the lower arm provided in our case by the interaction with the cantilever, while there is due to a gate 
voltage. 

In conclusion, we have analyzed the transmittance and the phase shifts behavior in a two-terminal AB ring whose 
lower arm is coupled to a cantilever. We showed that the phase rigidity can be broken as an effect of the inelastic 
scattering channels introduced by the cantilever, and a continuous phase shift can be obtained through the measure- 
ments of the linear conductance. The continuous phase variation as a function of the incident electron energy in 
experiments can be exploited to obtain the value of the electron-phonon coupling. Our proposal is within reach with 
today's technology employed in nanoelectromechanical systems. It can be realized by means of a semiconducting ring 
(e. g. In As^S) operating at mK temperatures with radius R ~ 60 4- lOOnm coupled to a molecular cantilever (e. g. 
a cantilevered or bridged single walled carbon nanotube^i) with a mass w 10~ 24 -4 10~ 23 kg and with fundamental 
frequency u> ~ 10 2 -j- 10 3 GH z, thus allowing to realize dwell times for the electrons less or of the order of the cantilever 
frequency. The cantilever length can be taken of the order of the (0.02-0.05)/im and its distance from the arm of the 
ring of the order (0.01-0.1)/im. In this way an extra electron in correspondence of the position of the tip causes a small 
vertical displacement. The phase rigidity breaking can be regarded a common feature of two-terminal nanomcchanical 
systems and thus we propose measurements of phase shifts as a way to determine the strength of the electron-phonon 
coupling in NEMS. 
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FIG. 1: Device: A mesoscopic ring symmetrically coupled to two external leads and a cantilever. 
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FIG. 2: Transmissions coefficients Tn (thick biack fine), T21 (full line), and %i (thick grey line) as a function of the Aharonov- 
Bohm flux / with the following choice of parameters: kL = 3.5, 7 = 0.2, pL = 3. The dashed curve represents the total 
transmission including also 3 and 4 phonons contribution. 
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FIG. 3: Transmissions coefficients 7n (thick black line), T21 (full line), and 7oi (thick grey line) as a function of kL for the 
following choice of parameters: 7 = 0.2, pL = 3 and / = 0.49, 0.5, 0.51 (top to bottom). 




FIG. 5: Phase shift 5i computed for: kL — 7 (box) and kL = 5.5 (diamond) [left panel] and kL = 8.9 (box) and kL — 9.9 
(diamond) [right panel] by fixing pL = 3. In all the curves 5i = tt when 7 = 0. 



